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T— I ' ABSTRACT 

^ ' We present a new set of stellar interior and synthesis models for predicting the inte- 

, grated emission from stellar populations in star clusters and galaxies of arbitrary age 

and metallicity. This work differs from existing spectral synthesis codes in a number of 
I important ways, namely (1) the incorporation of new stellar evolutionary tracks, with 

. sufficient resolution in mass to sample rapid stages of stellar evolution; (2) a phys- 

' ically consistent treatment of evolution in the HR diagram, including the approach 

(«~^ \ to the main sequence and the effects of mass loss on the giant and horizontal-branch 

. phases. Unlike several existing models, ours yield consistent ages when used to date 

' a coeval stellar population from a wide range of spectral features and colour indexes. 

We use Hipparcos data to support the validity of our new evolutionary tracks. We 
Q ■ rigorously discuss degeneracies in the age-metallicity plane and show that inclusion of 

|H I spectral features blueward of 4500 A suffices to break any remaining degeneracy and 

C/3 . that with moderate S/N spectra (10 per 20A resolution element) age and metallicity 

' are not degenerate. We also study sources of systematic errors in deriving the age 

^1 of a single stellar population and conclude that they are not larger than 10 — 15%. 

•j-^ . We illustrate the use of single stellar populations by predicting the colors of primor- 

' dial proto-galaxies and show that one can first find them and then deduce the form 

^ I of the IMF for the early generation of stars in the universe. Finally, we provide ac- 

■ - - I curate analytic fitting formulas for ultra fast computation of colors of single stellar 

populations. 

Key vifords: galaxies: stellar populations; stars: stellar evolution 



1 INTRODUCTION 

The synthetic stellar spectrum of a galaxy is a well 
established theoretical tool for investigating the proper- 
ties of the integrated light from distant galaxies where 
individual stars cannot be resolved. Since the late 
60's several groups have developed different grids of 
synthetic stellar population models using a variety of 
stellar interio r tracks and o bs erved or theoretica l stella r 
spectra Ce.g. [Tinslevl ^ll968^: iBarbaro fc BerteUil Jl977 1 



Renzinil lll98ll): iBruzuall JlQSa): iBarbaro fc OlivJ ]l98e 1 



_Arimoto & Yoshii ('l987');'Guiderdoni & Rocca-Volmeranee' 
11987.'): .Bruzual fc Chariot (.199a'): .Worthev. 099; 
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11994) : iFioc fc Rocca-Volmerangel J1997. . 
19981) : iLee et all ^l2002^ '). Theldea is very 
simple: stars are born with a given initial mass function 
(IMF) and they evolve in time according to stellar evolu- 
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tion. At any time one can compute the integrated spectrum 
by summing up the individual spectra of the stars in the 
population at that instant. Two major approaches have 
been used to compute the integrated ligh t of a stellar 
population: the fuel con sumption theorem llRenzinilll98ll : 
iRenzini fc Buzzonilll983h and the isochr one technique, first 
developed by iB e^aro fc, Bertelll il977! ) and later used by 
ICharlot fc Bruzual. LlddW The fuel consumption theorem 
simply uses the fact that the contribution of stars in any 
post main-sequence evolutionary stage is proportional to 
the amount of nuclear fuel they burn at that stage, and 
approximates the post main-sequence evolution of the stars 
in the integrated population by the stellar track of the most 
massive star alive at the main sequence turn-off (MSTO). 
In contrast, the isochrone technique uses a continuous 
distribution of stellar masses, and hence tracks, to compute 
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the locus in luminosity and T^s at a given time for any 
mass. From this, a smooth isochrone can be computed. 
The fuel consumption theorem remains an elegant method 
of studying the fastest stages in the evolution of any 
population. On the other hand, with the new generation of 
fast computers, the isochrone technique is clearly the most 
accurate and precise method of computing the integrated 
light of any stellar population. 

Regardless of the computational technique used to cal- 
culate the integrated spectrum of a stellar population, the 
most important ingredient remains the stellar input: both 
stellar interior and stellar atmospheric models. Disagree- 
ment over stellar interior models (convection, mass loss, 
opacities), the modelling of post main sequence evolutionary 
stages and the modelling of stellar atmospheres (opacities, 
NLTE effects, bolometric correction, mass loss, etc.), com- 
bine to make the derived age of even a simple (i.e., no dust, 
no AGN contamination) stellar population vary by about 
10% (for a fixed mass and metallicity), depending on which 
of t he currently available s ynthetic stellar popu lation codes 
fe.g. lCharlot et al.l Jl99d) : Ispinrad et al.l Jl997l) 'l is used to 
interpret the data (see section 4). A similar problem occurs 
when trying to date Globular Clusters in the G alaxy, the 
age of which is currently uncertain bv about 10% jchaboveil 
Il995l : ljimenez et al ] |l996l:IChabover fc Krausjl2002h . 

We have been motivated to attempt to improve this sit- 
uation by the fact that the new generation of large 8- 10m 
telescopes can now deliver spectra of galaxies at 2 > 1 of 
sufficient quality to merit accurate age dating. The accurate 
determination of the ages of high-redshift galaxies can yield 
important constraints not only on models of galaxy forma- 
tion, but also on t he age of the universe. In particular, in 
a series of papers JPunlop et al.lll996l : ISoinrad et al.lll997t 
iNolan et al.ll200lLl2003l) ' we have addressed the issue of de- 
termining the ages of the reddest known elliptical galaxies at 
2 ^ 1. To aid in the interpretation of our data and others, 
we have developed a set of simple synthetic stellar popu- 
lation models (SSP) that overcomes some of the problems 
described above. 

Most of the disagreement between the existing synthetic 
stellar evolution codes stems from the difficulty of modelling 
accurately the post main-sequence evolution, both because 
the physics involved in these stages is not completely known 
(opacities, convection, nuclear rates) and because mass loss 
strongly controls the final fate of the evolution of the star in 
these phases. Ideally, one would like to have a robust set of 
stellar models that are computed self-consistently (i.e. inte- 
rior, photosphere and chromosphere computed at the same 
time) and that include the effects of mass loss and dust grain 
formation. However this is not yet possible, and in any case 
it is important to realise that mass loss cannot be incorpo- 
rated as a fixed quantity for all stars in the population since 
it varies from star to star. 

In the new synthetic stellar population models pre- 
sented here we have endeavoured to improve the ability 
of the modelling to accurately reproduce the post-main se- 
quence evolution of real stellar populations by incorporating 
an algorithm which has been previously applied with suc- 
cess to a numbe r of other stellar evolution problems (e.g. 
I.Iorg-ensenI ri99l:I.Ior<xensen Thenlll ri99l:I.Timenez etltl 
l)l99ff ): Ijorgensen fc Jimened (|l223))- This algorithm accu- 
rately simulates the evolution of all post main-sequence evo- 



lution stages, and includes a proper modelling of the HB 
along with an accurate account of the formation of carbon 
stars on the AGB. Furthermore, the mixing length parame- 
ter and the mass loss are properly calibrated using the posi- 
tion of the RGB and the morphology of the HB in real star 
clusters, respectively. 

The other main new feature of our spectral synthe- 
sis models is the incorporation of new stellar evolutionary 
tracks, with sufficient resolution in mass to sample rapid 
stages of stellar evolution. As a result of these improvements 
(which we describe in detail in this paper), unlike several 
existing population codes, our models yield consistent ages 
when used to date a coeval population from a wide range of 
spectral features and colour indexes. 

Our new SSP code has been applied successfully to a 
variety of different populations, it has been used to deter- 
mine the ages of high redshift galaxies ( DunloD ot al. 199 ^, 
the ages of Low Surface B rightness Galaxies dPadoan et alJ 
ll997l : Tlimenez et al.l ll9 98[l. the role of star fo rmation and 
the TuUy-Fisher law (llIgaYgng_^^_jnrngngj ll999^ and the age 
of the Galactic disc i Jimenez et alJll998h . The purpose of 
this paper is to present the new library of synthetic stellar 
population spectra and discuss in more detail the physics 
and assumptions in our SSP modelling procedure. 

The paper is organised as follows: in §2 we present the 
new set of stellar interior tracks and discuss their accuracy 
when confronted with individual stellar observations. The 
synthesis models are presented in §3. The degeneracies in 
the age- metallicity plane are discussed in §4 while system- 
atic errors are considered in §5. The application of synthesis 
models to primordial proto-galaxies is presented in §6 along 
with a method to determine the initial mass function of these 
galaxies. §7 discusses the IR flux density and detectability 
of primordial proto-galaxies. Our conclusions are presented 
in §8. An appendix provides fitting formulas for computing 
broad band colors of SSPs. 

2 STELLAR MODELS AND PHYSICS INPUT 
2.1 Library of stellar interiors 

We have computed a new interior stellar library with the 
stellar evolution code JM STAR developed by one of us (JM) 
from the code of lEggleton (1971 , 1972.) . 

In this code, the whole star is evolved by a relaxation 
method without use of separate envelope calculations. Some 
of the advantages of this approach are that (1) gravother- 
mal energy generation terms are automatically included for 
the stellar envelope, (2) mass loss occurs at the stellar sur- 
face rather than an interior point of the star, and (3) the 
occurrence of convective dredge-up of elements produced by 
nucleosynthesis in the interior to the photosphere is clearly 
identifiable. The code uses an ada ptive mesh technique sim- 
ilar to that of lWerlev et alJ (Il984l) . Advection terms are ap- 
proximated by second-order upwind finite differences. Con- 
vective energy transport is modeled by using the mixing- 
length theory described by iMihalaai (1978). Composition 
equations for H, ^He, "^He, ^^C, "N,^'^0, and ^^Mg are 
solved simultaneously with the structure equations. Com- 
position changes due to convect ive mixing are modeled in 
the same way as lEggletonl il972h by adding diffusion terms 
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Figure 1. Stellar tracks in the Teff— L plane computed with JMSTAR are shown for 3 masses and metallicities. The tracks have been 
evolved from the initial Hayashi contracting phase until the white dwarf phase or carbon ignition, depending on the mass. The tracks 
are evolved in one run without interruptions and are able to go through the different He flashes in the AGB phase (see text for more 
details). 



to the composition equations. However, the prescription for 
the diffusion coefficient differs from that of EaEloton ( 1972). 
The di ffusion coefficient is consi stent with mixing-length 
theory dlben fc MacDonaldl Il995() . aeon = f3wconl, where 
uicon is the convective velocity, I is the mixing length, and 
/3 is a dimensionless convective mixing efficiency parame- 
ter. OPAL radiative opacities dlgiesias & Rogers 1996) are 
used for temperatures (in kelvins) above logjQ T = 3.84. 
For temperature below logm T — 3.78, we use opacities 



kindly supphed by D. Alexander and calculated by the 
method of Alexander & Ferguson ( 1994D • Between these 
temperature limits, we interpolate between the OPAL 
and Alexander opacities. Nuclear reaction rates are taken 
from lAn gulo et alJ (119 991) with screening corrections from 



ISalpeter fc van Horn' ('1969') and'ltoh et al] lll979D . Neutrino 
loss rates are from Bc audct ct al. ( 196 J) with modifications 
for neutral cu rrents (iRamad uraii Il976l) . Plasma neutrino 
rates are from iHaft et al .1 (Il994l) . The equation of state is 
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dete rmined by minimizat ion of a model free energy (see, 
e.g.. lFontaine et al] il977l) ') that includes contributions from 
internal states of the H2 molecule and all the ionization 
states of H, He, C, N, O, and Mr. Electron de generacy is 
included by the method of lEggleton et al.llll973l) . Coulomb 
and quantum corrections t o the equation of state follow the 
prescription of llben et alJ (Il992f) . with the Coulomb f ree en - 
ergy updated to use the results oflStrin gfellow et al.Ul99Cl) . 
Pressure ionization is included in a thermodynamically con- 
sistent manner by use of a hard-sphere free -energy term. 

Mass loss is included by using a scaled iReimer 
mass-loss law, 



M = r)Rl.27 X 10-''Kr^L' ''Tj'MQyr-^ 



(1) 



for cool stars (Tch ^ 10^ K) and an approximation to 
the theoretical result of lAbbofr3 il982l) . 



M 
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Meyr ' 



(2) 



for hot stars (Toft ^ 10'* K). For AGB stars we also 
include additional mass loss at a rate obtai ned by fitting 
the o bserved mass loss rates for Mira variables iKnapp et alJ 
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The surface boundary condition is treated as follows. 
We assume that the atmosphere is plane-parallel and thin. A 
small value of optical depth (typically r = 0.01) is assigned 
to the center of the outermost zone of the stellar model. 
The corresponding surface temperature Ts is related to the 
effective temperature, T^s, by the Eddington approximation 
Tf = rcff(0.5 + 0.75r). The surface gas pressure satisfies 
PgS — T(ges/k) where QeS is the effective gravity (surface 
gravity reduced by the effects of radiation pressure) at the 
surface of the star. We stress that the code treats the surface 
boundary on an equal footing with all other shells in the star. 

We have used the evolution code to create a grid of in- 
terior models with the following parameters. We computed 
stellar tracks for five different (solar scaled) metallicities: 
Z = 0.0002,0.004,0.01,0.02 and 0.05. We fixed the helium 
content for the solar model {Y — 0.27) as to reproduce 
the luminosity of the Sun at the present ag e and adopted 
dY/dZ = 2, consistent with the findings of l.Timenez et all 
l|200,^l . We show some representative stellar tracks in figQ 
for different metallicities. The tracks were evolved from an 
initial Hayashi phase to the very late stages in a continuous 
run. Thus it was possible to deal with the different flashes 
along the evolution during the giant phases. The evolution 
of low mass stars (M ^ SAf©) is ended after they settle onto 
the white dwarf cooling sequence. The evolution of massive 
stars (M ^ WMq) is stopped at the end of carbon burning 
or if the central temperature exceeds lO^-ft". For the interme- 
diate mass stars, the calculation is stopped once the envelope 
begins to be ejected hydrodynamically due to the radiative 
luminosity exceeding the local Eddington limit. 

In order to produce a model for the synthetic stellar 
population it is necessary to have a set of stellar interior 
tracks that includes all stellar evolutionary stages for a large 
range of masses. It is well known that the main discrepan- 
cies between different stellar population codes arise from 
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Figure 2. Evolution of a solar mass star with solar metallicity 
under the assumption that no mass loss occurs (left panel) or that 
it does {jj = 0.4; right panel). The dashed line is the red giant 
branch while the solid line is the asymptotic giant branch. The 
difference in luminosity, for the termination of the AGB, between 
the two cases is quite apparent (almost an order of magnitude). 
The diamond shows when a carbon star will be produced. By not 
including mass loss in stellar tracks one can seriously overestimate 
the luminosity of the giants in stellar populations. 
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0.0002 


0.235 


0.004 


0.24 


0.01 


0.25 


0.02 


0.27 


0.05 


0.33 



Table 1. Metallicity and helium fractions for the grid of stellar 
interior models computed with JMSTAR and used to assemble 
the grid of synthetic stellar population models. All stellar tracks 
were computed using solar-scaled abundances. 



the different treatme nt of the late stages of stellar evolution 
JCharlot et alll99dl . Since a proper and accurate modelling 
of the HB and AGB is extremely important for predicting 
the UV (HB) and IR (AGB) properties of any stellar popu- 
lation, we have included in our SSP code an accurate semi- 
empirical algorithm for computing the evolution of the HB 
and AGB. This is very important since evolution after the 
RGB cannot be described on the basis of one choice of mass 
loss parameter (this would result in the absence of a HB in 
globular clusters). The analytical description of the effect 
of mass loss (and the influence of different choices of mix- 
ing length) in our semi-empirical algorithm, makes it par- 
ticular useful in handling realistic population descriptions 
where variation from star to star in the mass loss efficiency 
needs to be described as a distribution function. We there- 
fore believe that we can obtain the highest accuracy by use 
of a semi-empirical algorithm. A direct comparison between 
the semi-empirical method and tracks computed with JM- 
STAR can be found in .Jimen ez ct al. (1995). The fast and 
accurate computation of the HB, AGB and TP-AGB, al- 
lowed by the semi-empirical method makes it possible to 
produce stellar tracks for several values of the mass loss 
parameter (rj in the Reimers formulation, or any other sim- 
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Figure 3. The morphology of the HB is shown for 4 different ages 
(all for solar metallicity) . We have use a gaussian distribution for 
the mass loss parameter {rj) around a central value of 0.4. Without 
this the HB would be a point. Note that for ages less than 10 Gyr 
the HB is actually vertical. 



ilar) and mixing-length parameter (a). Therefore, it is pos- 
sible to analyse the impact of these (unknown) parameters 
in the final stellar track (see Fig.|2l. In two prev ious studies 
Jjoreensen fc TheiUllTooa: I Jimenez et alJll99d) we demon- 
strated that the morphology of the HB is due to star-to-star 
variations in the mass loss parameter. This automatically 
gives a powerful tool to model the HB including, a priori, a 
realistic distribution of values for the mass loss parameter 
in the RGB. 

Mass loss takes place at the very end of the RGB and 
is most often described by the empirical Reimers mass loss 
law eimersl which is derived from measurements of 

M type giants. However, instead of following the evolution 
along the RGB using a fixed value for the rj parameter (as is 
usually done in the literature, and as was done by Reimers 
himself) we proceed in a different way. We use the computed 
grid of stellar interior models and calculate interpolation for- 
mulas (for the RGB) for L, T^b, M, Mc (core mass), Z, Y 
a nd a (mixing length). The fitting formulas can be found 
in ljimenez et all lll996f) and for brevity we will not repeat 
them here. The average value of rj is determined by matching 
the observed HB mass distribution. A very fast and accu- 
rate numerical computation of the evolution along the RGB 
can then be performed by taking advantage of the fitting 
formulas. The addition of mass to the core during a given 
time step in the integration along the RGB is determined 
on the basis of the instantaneous luminosity, the known en- 
ergy generation rate, and the length of the time step. The 
mass of the core (Mc) at the end of the time step deter- 
mines the new value of L according to the fitting formulas. 
The total stellar mass at the end of each time step is calcu- 
lated as the mass at the beginning of the time step minus 
the mass loss rate times the length of the time step. The 
evolution of the synthetic RGB track is stopped when Mc 



reaches the value determined in the fitting formulas for the 
He-core flash. A more detailed study of this approach and 
its comparison with other stellar models can be found in 
IJimenez etm lll995ll . 

Evolution along the AGB is much faster than the cor- 
responding evolution in the RGB and HB, lasting around 
5 X 10*' years. It has two differentiated phases: the early 
AGB (E-AGB) and the TP-AGB. The E-AGB occurs af- 
ter the helium core is completely converted to oxygen (and 
smaller amounts of carbon) and starts when the star is burn- 
ing helium in a hydrogen-exhausted shell. This phase ends 
when the hydrogen shell has passed through a minimum in 
luminosity. The TP-AGB consists of a phase of double-shell 
burning with thermal-shell flashes and strong mass loss. The 
modelling of the AGB evolution is rather complicated due 
to the above considera tion s . We have modelled it using the 
prescriptions in i JorgensenI lll99 J) and will not repeat them 
here. The most important issue is that we can compute when 
(or not) the star along its AGB will become carbon rich as 
well as whether the star will end the AGB as a planetary 
nebula or not. The advantage of this approach to modelling 
the AGB is that, as for the HB, we can produce accurate pre- 
dictions on the fate of the star and explore a large range of 
parameters for the mass loss, mixing length, chemical com- 
position, etc, in a fast and accurate way. The interpolation 
formulas are based on nu merical AGB com putations from 
JMSTAR, as described in I JorgensenI (Il99lll . and the evolu- 
tion is terminated when the core mass reaches the total mass 
of the star, or the total mass minus the mass of a planetary 
nebula. 

To summarize, the most important advantages of our 
method for isochrone synthesis are as follows: 

• A state-of-the art code to generate stellar interiors that 
can evolve the stellar track from the initial contracting 
gaseous phase to the flnal stages without interruption and 
that includes the latest advances in reaction rates and opac- 
ities. 

• A realistic distribution function for the mass loss effi- 
ciency parameter rj is determined on the basis of the HB 
morphology. The actual values of used in the present 
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work are taken from our previous studies of globular clusters 
Ijimenez et alJil99(J) . In particular, this distribution ensures 
a realistic integrated value of the blue and UV radiation from 
the SSP. 

• The mixing-length parameter is determined from fits to 
globular clusters with known metallicities. This is very im- 
portant for a realistic description of the red and IR radiation 
from the SSP models. 

• The semi-empirical algorithm has the same accuracy 
as the original numerical computations on which it is based 
(plus being improved in a and r] based on observations) , but 
it is computationally much faster. 

In order to illustrate the effect that changes in mass loss 
can have, we have computed, using the above prescription, 
the evolution of the Sun {Z = 0.02, Y = 0.28) with and with- 
out mass loss (using the Reimers formalism with r) = 0.4). 
In the first case the star evolves well into the AGB and de- 
velops into a carbon star before it will eventually transform 
through the post- AGB phase into a WD. In the second case 
the star leaves the AGB at a much lower luminosity (10 
times lower) and at a greater T^s (3000K instead of 2500K). 
The consequence for synthetic stellar population models are 
quite obvious: if mass loss is ignored the red-IR luminosity 
of the population is overestimated (see Fig. 

Stellar isochrones are computed by interpolating stellar 
tracks at equivalent s tellar evolutionary points as defined by 
ISchaller et al.l il992f) . This is important in order to have a 
good sampling of all stellar evolutionary stages. Some ex- 
amples of stellar isochrones are shown in Fig. |11 where for 
clarity we are not plotting the post-AGB phase. 

2.2 Library of stellar photospheres 

To compute the synthetic spectra corresponding to our 
stellar evolutionary interior models, we have used a set 
of theoretical stellar photospheres for single stars, most 
of them come from the most recent Kurucz library 
( http : / /kurucz . harvard . edu/ grids . html ) . The accuracy 
of the Kurucz models has been the subject of recent studies. 
They do an excellen t job in repro ducing broad band colors 
of individual stars jBesseHllTgoi) stars. Some comparison 
has been done also between the Kurucz models and the, low 
resolution, spectra of about solar metallicity star ('Wilcots' 
Il994lh with equally good results. The models are known to 
not fare that well for wavelengths below ~ 2500 A or for very 
low temperatures. For temperatures cooler than AQGQK we 
have generated our own (LTE) models using the the version 
of the MARCS code developed by Uffe Jorgensen (private 
communication). The above library of atmospheres is used 
to compute spectra and magnitudes (using the appropriate 
filters) for the stellar interior isochrones. We choose theo- 
retical spectral libraries as for them one knows exactely the 
value of Teff, logg and the metallicity, thus one can match 
exactely them to the interior library. This is not the case for 
observed spectra, for which the above parameters are not 
exactely known, and worse, none of them have exactely the 
same value as the one needed for the interior models. We 
therefore much prefer to build theoretical models (contrast 
them with individual stars) and use them to predict the 
spectral energy distribution from the corresponding interior 
model. 
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Figure 5. Comparison of some of our isochrones with appropri- 
ate stars from the Hipparcos catalog. The isochrones are from 
the left [Fe/H]=0.50, 0.00 and 0.30 and shown by dotted, solid 
and dashed, respectively. Crosses refer to the stars with the same 
metallicity as the metal- weak isochrones, open circles are the 
stars matching with the solar isochrones and squares are the 
stars matching with the metal-rich isochrones. The metal- weak 
and solar isochrones fit the data very well, while the metal-rich 
isochrone provides a less satisfactory fit for Aly < 6. The ages of 
the isochrones are 10, 12 and 14 Gyr. 



2.3 Validity of our stellar interior models 

Here we show that our stellar models fit observations of indi- 
vidual stars. The first test concerns stars of different metal- 
licities in the solar neighborhood with accurate distances, 
while the second test benchmarks the models against star 
clusters in the Milky Way. 

We have used t he sample of sta rs from the Hipparcos 
catalog described in iKotoneva et al ] 12002) and for which 
accurate metallicities have been derived either from spec- 
troscopy or narrow band photometry. The sample consists 
of 213 stars with metallicities ranging from —1.3 < [Fe/H] 
< 0.3. All stars have accur ate photometry and p arallax dis- 
tances from Hipparcos (see lKotoneva et al.l (|2003))- In Fig.|5| 
we show the locus of a set of isochrones from our stellar 
library for different metallicities and ages (at such low lu- 
minosities the age makes no difference in the locus of the 
isochrone). The different symbols correspond to Hipparcos 
stars from the above mentioned catalog and with the appro- 
priate metallicity. Both solar and sub-solar models provide 
a good fit to the data, while the super solar model is not 
so good at higher luminosities but provi des a fair fit to the 
data for Mv < 6. iKotoneva et alj ll2002h discuss how other 
independent sets of isochrones fit the data and conclude that 
our isochrones are among the best fit to the Hipparcos data. 
Despite this, attention should be paid to super-solar models 
since they seem to be difficult to construct. 

Having demonstrated that the isochrones provide a 
good fit to dwarfs of different metallicities we explore how 
successful they are at fitting the giant branches. Since the 
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Figure 6. Observed colour magnitude diagram for three stellar clusters. The solid line corresponds to the best fitting isochrone from 
our models. For 47Tuc the broken solid line is a fit to the observations. Note the good fit to the overall shape of the color-magnitude 
diagram, both the MSTO and giant branches. See text for more details. 



age does play a role in the temperature of the giant branch, 
the test with Hipparcos local stars is not possible. In order to 
attempt this test we are therefore forced to use star clusters 
which have much more uncertain distances (and therefore 
less certain ages as well) than Hipparcos nearby stars but 
have a simpler star formation history: they are single bursts. 
However the colour of the giant branches is well determined 
and therefore provides a good test for the models, specially 
at low temperatures. So we do adopt distances published in 
the literature for these cluster and fit the best isochrone to 
the overall color-magnitude diagram. 

In Fig.|S]we show fits to the color- magnitude diagram of 
several globular and open clusters. In all cases the isochrones 
provide good fits to the overall shape of the color-magnitude 
diagram. Further, both the red giant branches and main 
sequence turn off' are well fitted by the isochrones. The fit 
is good for different colours: the red giant branch of M67 
is weU fitted for B - V and V - K as well. The fit to the 
globular cluster 47 Tuc has been done for an age of 11 Gyr 
and [Fe/H] = -0.76. For M76 we used [Fe/H] = -0.1 and 
an age of 5 Gyr. Finally, NGC188 has been fitted with an 
isochrone with [Fe/H] = —0.1 and 6 Gyr. 

We conclude that the stellar interior models provide 
good fits to available photometry of single stellar popula- 
tions and individual stars. We are most encouraged by the 



good fit the isochrones provide to Hipparcos dwarfs since 
these are the only stars with accurate distances, and there- 
fore absolute magnitude, and metallicities and their color 
magnitude diagram does not depend on the age of the star. 
It is also reassuring that the colours of theoretical red giant 
branches are able to reproduce those of observed clusters. 



3 THE LIBRARY OF SYNTHETIC STELLAR 
POPULATION SPECTRA 

We define a stellar population as a collection of stars formed 
at a certain rate in a certain chemical environment and 
whose integrated light we are able to measure. We define 
a single stellar population (SSP) as a stellar population 
formed in a uniform chemical environment during a burst of 
infinitesimal duration (delta function). The steps involved 
in building a SSP are: 

(i) A set of stellar tracks of different masses and of the 
metallicity of the SSP is selected from our library. 

(ii) An isochrone is built for the corresponding age. Each 
point in the isochrone heis a luminosity, effective tempera- 
ture and gravity that are used to assign the corresponding 
photospheric model. 
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Figure 7. Spectral energy distributions of single stellar populations for three different metallicities: 0.2, 1 and twice solar (from left to 
right) and ages between lO'^ and 13 X 10^ years (from top to bottom in each panel). 



(iii) The fluxes of all stars are summed up, with weights 
proportional to the stellar initial mass function (IMF). 

SSPs are the building blocks of any arbitrarily compli- 
cated population since the latter can be computed as a sum 
of SSPs, once the star formation rate is provided. In other 
words, the luminosity of a stellar population of age to (since 
the beginning of star formation) can be written as: 



L\^^{Z, to 



t)dZdt 



(4) 



where the luminosity of the SSP is: 



Lf ^^(Z, to-t)^ / SF{Z, M, t)lx{Z, M, to - t)dM (5) 

and 1\{Z, M, to—t) is the luminosity of a star of mass M, 
metallicity Z and age to—t, Zi and Zf are the initial and final 
metallicities, Md and Mu are the smallest and largest stellar 
mass in the population and SF{Z, M, t) is the star formation 
rate at the time t when the SSP is formed. In this paper we 
concentrate in building SSPs, while in previous papers we 
discuss ste llar populations that are c onsistent with chemical 
evolution djimenez et al.lll998ll200ol) . 

We have computed a large number of SSPs, with ages 
between 1 x 10® and 14 x 10® years, and metallicities from 
Z = 0.0002 toZ = 0.1, that is from 0.01 to SZq. Fig.Ejshows 
a set of synthetic spectra of SSPs for three metallicities and 
ages between 1 x 10^ and 14x lO'' years (from top to bottom). 



4 PARAMETER DEGENERACIES: AGE AND 
METALLICITY 

SSP are simple, given an IMF, they only depend on two 
parameters: age and metallicity. It is therefore customary 
to ask the question: can one recover these two parameters 
from an observed spectrum? Of course, no single galaxy in 
the universe is a SSP, since its star formation history is not 
going to be a single burst, nor its metallicity a single one. 
However, most (if not all) elliptical galaxies have formed 
their stellar populati ons at high redshifts in a very short 
duration of time (e.g.. |jimenez et alJ (^9^) and they have 
relatively uniform chemical compositions, thus their spectra 
can be correctly described by a SSP. 

The so called age-metallicity degeneracy has been the 
focus of much attention during the past few years (e.g. 
IWilcotsI lll994l) ). Several authors have argued that broad- 
band colors are degenerate and therefore cannot be used 
to determine simultaneously age and metallicity of a SSP. 
Iwilcots (1994) has pointed out that certain absorption fea- 
tures can be used to better determine age and metallicity 
simultaneously. These feature are usually confined to a rel- 
atively narrow spectral range (4000 - 7000 A). 

In principle there is no need to discard information 
across the spectrum since all of it may contain informa- 
tion that yields a joint determination of the two parameters. 
Here we show how the whole spectrum contains information 
about both age and metallicity and that inclusion of light 
bluewards of 4500 A helps to break any degeneracy in the 
age-metallicity plane. 

We first select a spectrum from the grid for a given 
age and metallicity and add random noise. We then try 
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Figure 8. Joint likelihood contours for age and metallicity. The 
four panels correspond to spectra whose actual age is 1,4,8 and 
12 Gyr and solar metallicity with spectral range from 4500 to 
7000 A and with S/N = 10 per 20A . The contours show the 
68, 95 and 99% confidence levels of recovering the true age and 
metallicity. Note that there is some age-metallicity degeneracy, 
although it is more pronounced for ages older than 8 Gyr. 



1,10 



Jimenez (true age/Gyr) 1 4 8 12 
BC (recovered age/Gyr) 1.2 3.5 7 11.1 



Table 2. The age for a fiducial model created from our set of 
models by adding random noise with S/N = 10 per resolution 
element of 20A as recovered by the Bruzual &; Chariot models, 
this gives an estimate of how important systematics are in de- 
termining the ages of stellar populations. From the above table 
it transpires that systematics in the absolute age are only at the 
10% level, while for the relative age are only at the 2 — 3 % level. 

seen that age and metallicity are not completely degenerate 
since the contours do close, however the joint determina- 
tion of age and metallicity is accurate only at the 20% level. 
This also shows that one does not need to concentrate in 
special spectral features chosen empirically, the entire spec- 
trum (despite the limited spectral range considered) does 
contain information. 

The above mild degeneracy is completely removed if we 
add 2000 A of bluer spectral light. Fig.|3]shows similar plots 
as the previous figure, but in this case the spectral range of 
the spectrum goes from 2500 to 7000 A similar to that of 
the SDSS survey. The figure shows that in this case age and 
metallicity can be recovered with accuracy better than 1% 
(des pite the moderate S /N (10) of the mock test spectra). 

iDunlop et alJ i2003l) reach similar conclusions using real 
spectra. They find that addition of the same amount of UV 
light is sufficient to completely determine the age-metallicity 
and star formation history of an elliptical galaxy. 
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Figure 9. Same as fig. ISl but now the spectral range has been 
increased to 2500 — 7000 A . Note the total removal of the age- 
metallicity degeneracy. 



to recover the age and metallicity using x minimization, 
marginalizing over the amplitude. In the first case we limit 
the spectral range from 4500 to 7000 A and add random 
noise with S/N — 10 per resolution element (20 A). We do 
this for four ages (1, 4, 8 and 12 Gyr) and solar metallicity. 
The corresponding contour plots of the likelihood surfaces 
for 1,2 and 3a levels are shown in fig. |H1 Already, it can be 



5 SYSTEMATIC ERRORS IN SSPS 

In the previous section we have shown that given moderate 
5'/A'^ and wavelength spectral coverage, the parameters of 
a SSP can be accurately recovered. What is the impact of 
systematics in this parameter determination? 

To answer this question and also to assess how our mod- 
els compare to previous work, we have done the following: 
first we have compared o ur models to those by B ruzual 

6 Chario t (2000 mode ls; iBruzual fc Chariot! lll993f) 1 and 
Worthey iWilcot j|l994l) . Fig. 1101 shows the spectral energy 
distributions as a function of wavelength for 3 different ages. 
The 3 models differ among each other over the whole spec- 
tral range. The models are more similar at older ages than 
at ages of 1 Gyr. The typical rms variation between our 
models and those of BC is of 5 — 7% in intensity. How im- 
portant is this difference? To address this question we have 
performed the following test: we pick a spectrum from our 
models (with a fixed age and metallicity) developed in this 
paper and add gaussian noise with S/N = 10 per resolu- 
tion element of 20 A. We then attempt to recover the age 
using Bruzual & Chariot models. Results are presented in 
Table 0] As can be seen the systematic error in the absolute 
age is about 10% while the relative age difference is only of 
a 2 — 3%. We therefore conclude that although systematic 
differences remain between different models, their impact on 
parameter recovery is at the same level as the difference in 
the intensity of the models. 

We can also investigate the impact that mass loss can 
have in, for example, the colours of integrated single stellar 
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Figure 10. Comparison between different single stellar population models. The models are those by Bruzual & Chariot (BCOO), Worthey 
(W) and ours (SPEED). The three models differ among each other, albeit only at the 5 — 7% level in flux. The discrepancy at wavelengths 
smaller than 0.2 /im between the models are due to different treatment of post-AGB phases (which arc not included in W94 models). 



populations. In Table|31we show the change in the integrated 
colours that different prescritions for mass loss, which trans- 
lates mostly in the morphology of the HB and termination 
of the AGB, produce. The impact can be significant in the 
colors if mass loss is ignored (as much as 0.3 magnitudes in 
the infrared), thus careful attention must be payed to this 
source of systematics. 



6 SSP AT ULTRA-HIGH REDSHIFT: 
PRIMORDIAL PROTO-GALAXIES 

At 2: > 5, traditional optical bands {U, B ,V, R) fall below 
the rest frame wavelength that corresponds to the Lyman 
break spectral feature (1216 A), where most of the stellar 
radiation is extinguished either by interstellar or intergalac- 
tic hydrogen. Because of this, galaxies at z > 5 are practi- 
cally invisible at those photometric bands, and even if they 
were detected, their colors would provide very little infor- 
mation about their stellar population. New detectors and 



space telescopes, such as SIRTF (http://sirtf.caltech.edu ) , 
and in particular the JWST (http://www.stsci.edu/jwst), 
will offer the possibility of detecting very distant galaxies 
at IR wavelengths, and of using photometric rodshifts also 
with f a r-IR colo rs, as proposed by iSimpsori fc Eisonhard3 
^1999^ : lPanagi^et al] ll20o3) . 

Nearby star forming galaxies are known to contain a 
considerable amount of dust, that extinguishes a large frac- 
tion of their stellar UV light and boosts, even by orders of 
magnitude, their IR luminosity. The effect of dust must be 
taken into account in the computation of colors involving 
photometric bands that span a large wavelength interval. 
However, the search for very young proto-galaxies, perhaps 
the first significant star forming systems in the Universe, 
might be pursued without having to model the effects of 
dust. The first stars to form in proto-galaxies must have 
primordia l chemical composition , or at least very low m etal- 
licity (e.g. lTegmark et all Jl997ft : lPadoan et al.Ul997M . Al- 
though it is difficult to model the formation and disruption 
of dust grains on a galactic scale, it is likely that the dust 
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AiU - B) 


A(B - V) 


A(V-/) 


A{B-K) 


0.5 Gyr 


N-AGB 


-0.10 


-0.08 


-0.10 


-0.25 


R-AGB 


-1-0.10 


+0.08 


+0.12 


+0.30 


N-HB 


-0.01 


-0.05 


-0.10 


-0.20 


B-HB 


-0.05 


-0.02 


-0.05 


-0.05 


R-HB 


-1-0.05 


+0.02 


+0.05 


+0.04 


5 Gyr 


N-AGB 


0.00 


0.00 


-0.05 


-0.10 


R-AGB 


0.00 


0.00 


+0.05 


+0.20 


N-HB 


-1-0.01 


0.00 


-0.04 


-0.20 


B-HB 


-0.03 


-0.02 


-0.03 


-0.05 


R-HB 


0.00 


0.00 


+0.03 


+0.03 


14 Gyr 


N-AGB 


-1-0.01 


+0.01 


+0.01 


-0.08 


R-AGB 


-1-0.05 


+0.03 


+0.10 


+0.20 


N-HB 


-0.10 


-0.04 


+0.04 


+0.10 


B-HB 


-0.15 


-0.05 


-0.05 


-0.2 


R-HB 


-1-0.05 


-0.05 


-0.05 


-0.05 


Table 3. The influence 


on the computed colours for 


a single syn 



thetic stellar population of changes in the assumptions underlying 
the population computations. We list the changes corresponding 
to a reference stellar population computed witha Salpeter IMF 
and solar metallicity. N-AGB corresponds to models computed 
without an AGB, while R-AGB corresponds to models computed 
with AGB but no mass loss. N-HB stands for models computed 
with no HB. R-HB and B-HB are models computed with red 
[r] = 0.0) and blue (t; = 0.7) HBs. The different approaches can 
change the predicted colours by up to 0.3 magnitudes. 

content of a galaxy grows together with its metallicity, and 
that a proto-galaxy with primordial chemical composition 
or very low metallici ty (Z ^ O.OIZ©) has practically no dust 
(e.g. Jimenez ct al] (|l99fll )'l. In this work we refer to young 
protogalaxies with metallicity Z ^ 0.01^0 as primordial 
protogalaxies, or PPGs. We show that, because of their very 
blue colour, PPGs do not suffer from very strong cosmolog- 
ical dimming, a feature which should aid their detectability 
in deep infrared surveys with the JWST. 

6.1 Spectral energy distribution and IR colors of 
PPGs 

We have computed the SED of PPGs using the extensive 
set of synthetic stellar population models developed in this 
paper. In particular, PPGs are modelled using a constant 
star formation rate and Z — 0.01.^0 metallicity. As noted 
above, we expect (and therefore assume) that no photomet- 
rically significant quantity of dust has already been formed 
in PPGs. We also assume that PPGs are young objects (age 
less or equal than 100 million years). Furthermore, and in or- 
der to test the shape of the primordial IMF, we have adopted 
a Salpeter IMF {x — —1.35) with four different low mass 
cutoff values: 0.1, 2, 5 and 10 M0 (i.e. the IMF does not 
contain any stars with masses below the cut-off values) . 

The effect of the different low mass cutoff values on the 
SED of a PPG is shown in fig. llll where the flux density F^ 
is plotted in arbitrary units for a 10 million year old PPG 
with Z — O.OIZq and constant star formation (solid lines). 
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Figure 11. Spectral energy distribution for a starburst with 
constant star formation rate, age of 10 million years, metallicity 
Z = O.OIZ0 and a Salpeter {x = -1.35) IMF with four different 
lower mass cutoffs: 0.1, 2.0, 5.0 and 10.0 M0 (solid lines), from 
top to bottom. Starbursts with no low mass stars (< 5 M0) 
are about a factor of 2 fainter than starbursts with low mass 
stars (> 0.1 M0), at 3 fim. Note that for cutoff masses larger 
than 5 M0, the spectral energy distribution changes very little, 
making our predictions almost independent of the exact value of 
the cutoff mass. The models have been normalis ed at 1000 A. 
Also s hown is a solar metallicity starburst from ISchmitt et alj 
Il997|) (dashed line). We finally plot an idealized starburst model 
with constant star formation rate, Salpeter IMF, age of 10 million 
years, Z = Q.2Zq and no dust (dotted line). The UV SED of 
this starburst can be a rough approximation of the UV SED of 
quasars. 



Different solid lines corresponds to the different IMF cutoffs, 
0.1, 2, 5 and 10 Mq, from top to bottom (the 5 and 10 Mq 
are almost overlapped). As expected, the lack of low mass 
stars translates into a deficit of fiux at IR wavelengths. Fur- 
thermore, for a low mass cutoff larger than 5 M0 the SEDs 
do not differ much since stars with masses larger than 5 Mq 
have similar spectral energy distributions in the IR. The 
difference in relative fiux density over the wide wavelength 
range illustrated in fig. 1111 provides the means to test the 
IMF in PPGs. The SED of an observed nearby starburst, 
from lSchmitt et'al] lll997l) . is also plotted in fig. 1111 (dashed 
line), together with a 10 million year model starburst, with 
continuous star formation, metallicity Z = 0.2Zq, and no 
dust extinction or emission (dotted line). The difference in 
the spectral slope between PPGs and nearby starbursts is 
striking, and it is still significant between PPGs and the 
idealized starburst model with no dust. Although the ideal- 
ized starburst model with no dust is unlikely to describe any 
galaxy observed nearby, it can be used as an approximate 
model for the UV SED of quasars, due to its rather fiat SED 
in the far ultra violet wavelengths. 

We have computed AB magnitudes (ttiab = 
— 2.51og(F^) — 48.59, where is expressed in erg s~^ Hz^^) 
for 3 different wavelengths: 1.2, 3.6 and 8 /im. Fig. ll2l shows 



12 Jimenez et al. 



the colour-colour trajectory for PPGs with a Salpeter IMF 
with a low mass cutoff of 0.1 Mq (thin solid line), and 5 
M© (thick solid line), metallicity Z — O.OlZe, and an age 
of 10 million years. The dashed line is the evolution of the 
nearby starburst from fig. 1121 and the dotted line the evo- 
lution of the idealized starburst model with no dust and 
Z = 0.2Zq, also from fig. 1121 All models have been sim- 
ply K-corrected and the numbers that label the trajecto- 
ries indicate the corresponding redshift. The figure shows 
that PPGs with redshift range 5 < z < 10 have colors in 
the range of values —3.0 ^ (1.2/im— 8/im)AB ^ —2.0 and 
—0.9 ^ (3.6/im— 8/im)AB ^ —1.5, that are inside the grey 
area. The idealized starburst model with no dust and sub- 
solar metallicity enters marginally the dashed area, but only 
for low redshifts {z < 0.5). This idealized model is unlikely 
to describe nearby galaxies, while it can be a rough descrip- 
tion of the UV SED of quasars, in which case it could be 
concluded that PPGs should have very different colors than 
quasars at redshift z > 0.5 (at least about 0.5 mag away 
in the color-color plot of fig. I12II . Nearby galaxies are ei- 
ther star forming galaxies with significant amount of dust 
(irregular, spiral or starburst galaxies), or older stellar sys- 
tems with little gas or dust (elliptical galaxies). In both cases 
nearby galaxies are much redder than PPGs. In fig. 1121 the 
dashed dotted lines show the color-color redshift evolution 
of two typical nearby galaxies (a spiral and an elliptical 
galaxies, from ISchmitt et al.l ijl^Ssi)). These galaxies, even 
if nearby, are always at least 2 mag redder than PPGs in 
the (1.2/xm— 8/xm)AB color. 

The fact that PPGs are the bluest objects, in this 
IR color-color plot, is an important result, since a photo- 
metric search for high redshift galaxies would, in princi- 
ple, be biased toward selecting the solar metallicity star- 
bursts, which arc the reddest galaxies, as already proposed 
by|SimDSon & Eiscnhardt (1999). Although it cannot be ex- 
cluded that galaxies of metallicity close to the solar value 
might exist at z = 10, and that PPGs are rare (since they 
are young by definition), fig. 1121 shows that primordial star 
formation at high redshift should be searched for in very 
blue objects. The shaded area in fig. I12l marks the expected 
location of PPGs. 



7 THE IR FLUX DENSITY OF PPGS 

We have just shown that PPGs can be photometrically se- 
lected as the bluest galaxies in the Universe. The question 
that needs to be answered now is: Can PPGs be detected 
at all with future telescopes such as the SIRTF and the 
JWST? PPGs could in fact be very faint because they could 
be very small, or because their star formation rate (SFR) 
could be very low. To answer this question, we have com- 
puted the expected fiux in njy per unit of SFR (in Mq/jt), 
as a function of redshift, for a PPG with a Salpeter IMF and 
cutoff mass of 5 Mq, in a A dominated Universe {^Im ~ 0.3, 
A = 0.7, Ho = 65 km s"^ Mpc"^). The plots are roughly 
independent of the starburst age, for age larger than a few 
million years. Fig. ll3l shows that, although PPGs do not ex- 
hibit a negative K-correction as galaxies do in sub-mm and 
mm bands, they suffer from very little cosmological dim- 
ming, even in an open Universe. From fig. 1131 it can be seen 
that PPGs with SFR of about 100 M© /yr have a flux in the 
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Figure 12. Colour-colour (in the AB system) redshift evolu- 
tion of; the observed starburst of solar metallicity from Figure 1 1 
(dashed line); the idealized starburst model of 1/5 solar metal- 
licity and no dust, also from Figure 11 (dotted line); two PPGs 
models (from Figure 11), with constant star formation rate, age 
of 0.01 Gyr, Z = O.OIZq, and Salpeter IMF with 0.1 Mq cutoff 
(thin solid line) and 5 Mq cutoff (thick solid lin e); two typical 
nearby galaxies (a spiral and an elliptical), from ISchmitt et alj 
(1997). The shaded area shows where PPGs are expected to be. 



faintest band (8/xm) of about 1 nJy (although they would 
be much brighter at 1.2 /xm - about 10-20 nJy). The IRAC 
camera at SIRTF will be able to measure in 1 hour expo- 
sure a flux of 500nJy at the 5a level, and will not be able 
to detect PPGs which a SFR of less than 100 Mq yr"^ but 
much higher. On the other hand, the JWST time calculator 
(htpp:/www. stsci.edu/jwst) gives (for a point source using 
the MIR-ACCUM instrument with a resolution of 3, S/N=5) 
that exposing 300 hours one can get a fiux of about 5nJy at 
8^m. In order to distinguish PPGs we also need observations 
at 3.6 and 1.2 (see fig. ll2ll . using the same calculator we 
obtain (this time NIR-ACCUM, resolution 3, S/N=5) that 
it takes 3.63 hours to achieve 1 nJy at 3.6 fim, while it takes 
3.12 hours at 1.2 ^m. It is worth noting though that PPGs 
are expected to be about 1 order of magnitude brighter at 
1.2 and 3.6 than at 8 /im. Therefore, for the 4 nJy constraint 
imposed by the 8 /im band one would need only a few min- 
utes to achieve the required sensitivity at 1.2 and 3.6 /im. 
It is therefore accurate to say that the limit on our obser- 
vations comes from the 8 fim band solely. Now, inspection 
of fig. 1131 shows that a fiux of 5nJy at 8^m is achieved for 
PPGs with SFRs of about 400 Mq yr"^ For example, one 
needs to have systems with ages of 10^ yr and masses of a 
few 10^ Mq in order for them to achieve a fiux of 4nJy. This 
is not an unreasonable mass for the first objects since it is 
likely that star formation of the first o bjects will be retarde d 
until a relatively large halo is in place jOh fc Haima J2003l) . 
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stellar win ds that would stop the mass accretion onto the 
proto-star iAdams fc Fatuzzol[l993) : iii) fractal mass distri- 
bution with fragmentation down to one Jeans' mass fLarsorJ 
il992) . If gravitational su b-frag mentation during collapse 
is not very efficient (see iBosj l|l993)), the value of the 
Jeans' mass d e termi nes the lower mass cutoff of the IMF. In 
IPadoan et al.l il997l) , numerical simulations of super-sonic 
and super- Alfvenic dPadoan fc Nordlundl Il999l) magneto- 
hydrodynamic turbulence are used to compute the prob- 
ability density function of the gas density, which is used 
to predict the distribution of the Jeans' mass in turbu- 
lent gas, under the reasonable assumption of uniform ki- 
ne tic temperature. The Jeans' mass distribution computed 
in IPadoan et alJ il997l) has an exponential cutoff below a 
certain mass value, that is found to be: 




Figure 13. The expected IR flux density of PPGs, per unit of 
star formation rate, as a function of redshift, for four wavelengths 
and 2 cosmologies (in both cases Hq = 65 km s~l Mpc~^). Note 
that the cosmological dimming is very small since PPGs have 
Fu oc v^'^ (see Figure 11). The same model as in Figures 11 and 
12 has been used, with a 5 Mq cutoff IMF. The flux density has 
been computed for an age of 0.01 Gyr, and increases rather slowly 
afterward. 

7.1 The Primordial Stellar IMF 

Observational and theoretical arguments suggest that stars 
forming from gas of very low metallicity (< 10~^Zq) could 
have an initial mass function (IMF) shifted toward much 
larger masses than stars formed later on from chemically 
enriched gas. The observ ational a rgum ents are extensively 
discussed in a paper by iLarsorJ lll99^ . and we refer the 
reader to that work. The main theoretical reason in fa- 
vor of a 'massive' zero metallicity IMF is the relatively 
high temperature (T > 100 K) of gas of primordial com- 
position, that is cooled at the lo west te mperatures mainly 
by H2 molecules llPalla et aljll983 : Mac Low fc Shullll98a 
Shapiro fc Kane 1987: Kane ct^ Il99(i iKang fc Shapird 
1992. : .Anninos fc Norman 1996'). One of the most exciting 
aspects of the photometric discovery of PPGs would be the 
possibility of investigating the nature of their stellar IMF, 
once their redshifts are available. 

It is likely that while the local properties of the ISM do 
not interfere significantly with the self similar dynamics that 
originate the stellar IMF, they do play a role in setting the 
particular value of the mass scale where the self similarity 
is broken. According to this point of view, one expects the 
stellar IMF to have always more or less the same power law 
shape, down to a cutoff mass whose value depends on local 
properties of the ISM, and up to the largest stellar mass, 
whose value is limited either by the total mass of the star 
formation site, or by some physical process that prevents the 
formation of super-massive stars. 

The lower mass cutoff of the stellar IMF has been pre- 
dicted in models of i) opacity limited gravi t ation al fragmen- 
tation JSilklll977allbH:lYoshn fc Saidll985lll98l : ii) proto- 



where T is the gas temperature, n the gas density, and 
the gas velocity dispersion. Using the ISM scaling laws, ac- 
cording to which n^''^(Ji, ~ const, one obtains: 

M^i„«0.1Me(^)' (7) 

that is a few times smaller that the average Jeans mass (the 
Jeans mass corresponding to the average gas density), and 
therefore an important correction to more simple models of 
gravitational fragmentation, that do not take into account 
the effect of super-sonic turbulence on the gas density dis- 
tribution. 

If the ISM has primordial chemical composition, and 
the main coolant is molecular hydrogen, a temperature be- 
low 100 K is hardly reached, and the stellar IMF might 
have a lower mass cutoff of about 10 M m. S i milar l ower 
mass cutoffs arc obtained in the models bv lSilld Jl977di and 
Yoshii fc Saio (1986), who estimated typical stellar masses, 
based on molecular hydrogen cooling, of approximately 20 
and 10 Mq respectively. More recent numerical simulations 
of the collapse and cooling of cosmological density fluctua- 
tions of large amplitude (the first objects to collapse in the 
Universe), yield even larger values of the Jeans' mass, of 
the order of 100 M0 (Bromm, Coppi, fc Larson 1999; Abel, 
Bryan, fc Norman 1998) 

The discovery of PPGs could shed new light on the 
problem of the primordial IMF. The redshift evolution of 
the two colors (1.2/^m— 8/^m)As and {3.6fim~8fim) ab , com- 
puted with the PPG model discussed in this work, is plotted 
in fig. 1141 The solid line is the case of a PPG with a Salpeter 
IMF with 5 Mq cutoff, and the dashed line the same PPG 
model, but with a 0.1 Mq cutoff. Once a PPG candidate is 
selected with the IR broad band photometry as a very blue 
object (colors inside the shaded area in fig |12ll . and its red- 
shift is estimated with a Lyman drop method, or with an Ho? 
search, the IR colours provide a tool to discriminate between 
a standard IMF, and an IMF deprived of low mass stars. The 
lower panel of fig ll4l shows that a PPG with a 'massive' IMF 
can be about 0.5 mag bluer in (1.2^m— 8/^m)AB than a PPG 
with a standard IMF. 
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Figure 14. Colour evolution (AB system) with redsliift for PPGs 
with a Salpeter IMF with 0.1 Mq cutoff (dashed line), and 5 Mq 
cutoff (solid line). As expected from Figure 11, the biggest differ- 
ence occurs for colours that sample both the near-IR (rest— frame 
UV) and the far-IR. The difference between the Salpeter and the 
'massive' IMF is about 0.5 mag in (1.2/im— 8^m)^s . Once PPGs 
have been found at high z using the broad band color-color selec- 
tion from Figure 12, and their redshift has been determined with 
narrow band photometry, the (1.2/im— S/im)^^ color can be used 
to discriminate among objects with a standard Salpeter IMF or 
a 5 Mq cutoff IMF. 



needed to find all PPGs in the redshift range 5 ^ 2: ^ 10. 
Note that the equivalent width of Lyman-Q and Hq has 
been over-estimated here. Moreover, the advantage of deep 
broad band photometry is that, together with detecting and 
selecting PPGs, it provides at the same time important in- 
formation about their stellar populations. However, it is im- 
portant that the emission line technique (or a Lyman break 
technique) is available aboard the NGST, since photomet- 
ric redshifts measured with narrow filters will probably be 
the best (or the only) way to further constrain the redshift 
of PPG broad band photometric candidates, which is nec- 
essary to extract information about their stellar population 
from the broad band colors (fig. 4)^. 

A star formation rate of 100 M0/yr over a few million 
years is necessary for detecting a PPG with the JWST. With 
such SFR, after 10^ years 10^ Mq of gas is turned into stars. 
In order to still have a metallicity oi Z ^ O.OIZq , these stars 
must be formed in a system with baryonic mass of at least 
1 X 10^^ Mq, that is a very large galaxy, or a small group of 
galaxies. Such massive systems are inside large dark matter 
halos that are not collapsing yet at redshift 5 ^ « ^ 10. It 
is possible that PPGs that can be detected with the JWST 
are the progenitors of very large galaxies, in a phase when 
their dark matter halo has not turned around yet. If PPGs 
are discovered, their spectro-photometric properties could 
give very important clues for the problem of star formation 
in galaxies (such as the origin of the stellar IMF) and their 
luminosity, abundance, and redshift distribution would trace 
the complete history of the very first star formation sites in 
the Universe. 



7.2 Detectability via broad- band infrared imaging 
versus emission-line searches 

In this work we propose to select PPGs as the bluest ob- 
jects in deep IR surveys, on the basis of the color-color 
plot shown in fig. 2. We now address the question of how a 
broad band photometric selection of PPGs performs, com- 
pared with searches of emission lines, such as Lyman-Q and 
Hq. The rest frame equivalent widths of Lyman-Q and Hq 
can be very roughly estimated by assuming that each pho- 
ton below 1251 and 1025 A will originate a Lyman-Q and 
Hq photon respectively. The equivalent widths estimated in 
this way are of course upper limit to the true equivalent 
widths. We find that the equivalent width of Lyman-Q is 
380 A while the equivalent width of H— q is 4400 A -the 
latter is so high due to the fact that the continuum at 6563 
A is rather faint in PPGs (see Fig. IIH . Assuming that the 
lines have intrinsic widths at rest-frame typical of the virial 
velocity of a galaxy (for example a line width of 300 km/s 
corresponds to 2 and 13 A respectively), one finds that they 
will only be about 30 times brighter than the continuum at 
z ~ 10. Since one would need to shift the narrow filter for 
about 500 steps, or more, to search for all possible emitters 
in the redshift range 5 ^ z ^ 10, the advantage of the lines 
being brighter than the continuum is offset by the number 
of steps needed to find all PPGs between z = 5 and 10. It 
seems therefore that IR broad band photometry is an eas- 
ier way to both detect and select PPG candidates than the 
emission line technique, because only one deep exposure is 



8 CONCLUSIONS 

In this paper we have presented new stellar interior tracks 
and single stellar population models of arbitrary age, metal- 
licity and IMF. Our main findings are: 

(i) A new set of stellar interior models has been presented. 
The models cover all stages of stellar evolution and we have 
built isochrones out of them for ages between 10^^-1.4x10^°. 

(ii) It is possible to construct stellar evolution models 
that accurately reproduce the properties of individual stars 
for a wide range of ages and metallicities. 

(iii) We have presented a new algorithm to compute the 
evolution of stars in the RGB, HB and AGB. This algorithm 
makes it possible to explore the effect of variations in some 
of unknown parameters of stellar evolution like mass loss, 
mixing length etc. 

(iv) We have developed a new and fast algorithm to build 
synthetic stellar evolution spectra and colour-magnitude di- 
agrams of arbitrary metallicity and age. 

(v) We have shown that changes in the values of the stel- 
lar parameters like mass loss and mixing length can change 
the predicted colours of a population by as much as 0.4 mag. 



^ We note in passing that the acquisition of photometry around 
the Lyman break to determine their redshift would not be time 
consuming as shown before and that the biggest limitation in 
detecting PPGs is due to the poor performance of the 8 /^m de- 
tectors. 
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(vi) We have studied degeneracies in the parameter space 
(age and metalhcity) and shown that these parameters are 
only degenerated if the wavelength range of the spectrum is 
very small or only a few spectral features are chosen. Addi- 
tion of light bluewards of the 4500 A (2500 - 4500 A) sig- 
nificantly reduces this degeneracies and, in fact, lifts them. 

(vii) It has been shown that systematic errors among dif- 
ferent models are at the level of 10-20%, despite using com- 
pletely different stellar input physics. It should be possible 
to reduce these errors even further. 

(viii) We have studied the photometric properties of very 
young proto-galaxies with primordial or very low {Z = 
O.OIZq) metallicity and no significant effect of dust in their 
SED. We have named these galaxies "primordial protogalax- 
ies" , or PPGs. Using the methods of synthetic stellar popu- 
lations, we predict that PPGs are the bluest stellar systems 
in the Universe. They can therefore be selected in color- 
color diagrams obtained with deep broad band IR surveys, 
and can be detected with the JWST, if they have a SFR 
of at least 100 Mq/jt, over a few million years. We have 
discussed the possibility of using the IR colours of PPGs to 
constrain their stellar IMF, and investigated the possibility 
that the stellar IMF arising from gas of primordial chemical 
composition is more "massive" than the standard Salpeter 
IMF. Finally we have argued that broad band photometry 
can be more efficient than emission line searches, to detect 
and select PPGs. 

The models are available on the 
world wide web (www.roe.ac.uk/~jsd and 
www. physics. upenn.edu/'^raulj). We provide the stellar 
interior tracks presented in this paper and the single stellar 
population models and tools to compute photometry and 
synthetic stellar populations with arbitrary star formation 
histories. 
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9 APPENDIX: ANALYTIC FITS FOR 
COLOURS OF SSPS 

The magnitudes for a SSP (normalized to 1 Mq) as a func- 
tion of age and metallicity, for a given pliotometric band 
UBVRIJK, are approximated to within 4% by: 



Mx = - 

where 
X = 

Y = 
and 



-2.5x^^TO(i + l,j + l)F^ 



i=0 j=0 



5.76 + 3.18 log r + 1.26 log^ r + 2.64 log^ r 
+1.811og''r + 0.381og^ r, 



2.0 + 2.059 log C + 1.041 log^ C 
+0.172 log* C - 0.042 log"* C, 



(8) 



(9) 



(10) 



T = 



(11) 

from 



— c = — 

Gyr' Ze 

Luminosities are obtained simply 

{5.61, 5.48, 4.83, 4.34, 4.13, 3.72, 3.36, 3.30, 3.28} for 
{U, B, V, R, I, J, H, K, L}. The i and j values appear 
as exponents of X and Y, respectively, and as indexes 
defining elements of the Cx matrices, given by 
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